The Functional Renormalization Group and 0(4) scaling 
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The critical behavior of the chiral quark-meson model is studied within the Functional Renormal- 
ization Group (FRG). We derive the flow equation for the scale dependent thermodynamic potential 
at finite temperature and density in the presence of a symmetry-breaking external field. Within 
this scheme, the critical scaling behavior of the order parameter, its transverse and longitudinal 
susceptibilities as well as the correlation lengths near the chiral phase transition are computed. We 
focus on the scaling properties of these observables at non- vanishing external field when approaching 
the critical point from the symmetric as well as from the broken phase. We confront our numerical 
results with the Widom-Griffiths form of the magnetic equation of state, obtained by a systematic 
e-expansion of the scaling function. Our results for the critical exponents are consistent with those 
recently computed within Lattice Monte-Carlo studies of the 0(4) spin system. 

PACS numbers: 



I. INTRODUCTION 

Understanding the phase structure and the critical 
properties of strongly interacting matter is one of the 
central problems addressed in studies of QCD thermo- 
dynamics. Lattice Gauge Theory (LGT) calculations at 
finite temperature show a clear separation between the 
confincd_(hadronic) and deconfined (quark-gluon plasma) 
phases [J. 

Universality arguments imply that two- flavor QCD ex- 
hibits a second order chiral phase transition belonging to 
the same universality class as the 0(4) spin system in 
three dimensions 

H H EL IS i 0, 1 • Consequently, the 
long-range properties of the chiral phase transition can be 
explored within effective models of the same universality 
class, independent of the specific dynamics. 

One such model is the chiral quark-meson model, 
which often is used as an effective realization of the low- 
energy sector of QCD @. Thus, by studying the ther- 
modynamics of this model near the phase transition, one 
can explore the leading singularities of thermodynamic 
quantities in two-flavor QCD at the second order chiral 
phase transition. In other words, we can find the univer- 
sal properties of dynamical chiral symmetry breaking at 
finite temperature and/or density. However, at this point 
we stress that a description of the QCD phase transition 
within effective models is subject to limitations. First, 
in chiral effective models one cannot address the decon- 
fmcmcnt phenomenon in QCD and second, one cannot 
obtain quantitative information on the thermodynamics 
outside the chiral critical region. More generally, non- 
universal quantities depend on the specific dynamics of 
the theory and can therefore not be computed reliably 
within effective models. 

As a rule, phase transitions and critical phenomena 
are related with singular behavior of some susceptibility 
which in turn provide a measure of the fluctuations of 
a physical quantity. At a second order phase transition, 



there are long-range critical correlations that appear due 
to the presence of massless modes. The long-range corre- 
lations imply divergent susceptibilities and consequently 
large fluctuations. 

In the mean field approximation, which is often used 
to describe the chiral phase transition [n| EH, E2 , one 
fails to obtain a quantitatively correct description of crit- 
ical phenomena, since the influence of fluctuations and 
non-perturbative effects near the phase transition are ne- 
glected. However, methods based on the renormalization 
group (RG) can account for both of these important ef- 
fects (see e.g. (HI). The advantage of the RG method 
is that it can describe physics across different momen- 
tum scales. In particular, within the RG framework one 
can capture the dynamics of the long-range fluctuations 
near the critical point. Wilsonian RG techniques, such 
as the functional renormalization group (FRG). are par- 
ticularly useful in describing phase transitions [141 ] . One 
of the important results of RG theory is that it provides 
a rationale for universality, i.e., the fact that models can 
be grouped into universality classes that, depending only 
on the dimension of the system and the symmetry of the 
order parameter, show the same critical behavior. 

In this paper we use the FRG method in the quark- 
meson model to explore the universal properties of the 
chiral phase transition. This model has been studied us- 
ing different RG approaches @, GJ, EE E3, E3 both in 
vacuum as well as in medium, at finite temperature and 
chemical potential. Various critical exponents have also 
been computed [H, HE HE H3|- I n this paper we go 
beyond previous RG studies of this model. We analyze 
the critical region in the broken and in the symmetric 
phase close to the critical point and discuss the scaling 
functions for different physical quantities as well as the 
magnetic equation of state (23| . Furthermore, we com- 
pute the critical exponents for the order parameter and 
its transverse and longitudinal susceptibilities as well as 
for the correlation lengths. We also calculate the effec- 
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tive critical exponents [24[ in the presence of the chiral 
symmetry breaking field. Finally, we explore the scal- 
ing behavior of different ratios of susceptibilities and the 
order parameter [25| . 

We also show that within the FRG method, with a 
suitable truncation of the FRG flow equations, the quark- 
meson model exhibits a universal behavior that is in 
quantitative agreement with that recently found in lat- 
tice calculations of the 0(4) spin system in three dimen- 
sions [2^ |. 

The paper is organized as follows: In Sec. [II] we intro- 
duce the chiral quark-meson model and the FRG method. 
The flow equations for the scale dependent thermody- 
namic potential at finite temperature and chemical po- 
tential are also derived in this Section. In Sec. IIIII wc 
discuss the universal scaling properties of various phys- 
ical quantities and extract their critical exponents. In 
Sec. II VI we summarize our results and discuss their rele- 
vance. 




t 



FIG. 1: The order parameter as a function of the reduced 
temperature t for several values of the external field fi- 



ll. THE CHIRAL QUARK-MESON MODEL 
AND THE FRG FLOW EQUATION 

We employ the chiral quark-meson model to explore 
the scaling properties and the critical equation of state 
within the FRG method. This model is relevant for 
studying strongly interacting hot and dense matter with 
two degenerate light-quark flavors, since it is expected 
to belong to the same universality class as QCD at the 
chiral phase transition. 

The chiral quark-meson model represents an effective 
low-energy realization for spontaneous chiral symmetry 
breaking at the intermediate momentum scale 4^/^ ps I 
GeV. One can view this model as an effective model of 
QCD, where the gluon degrees of freedom have been in- 
tegrated out. Consequently, the model does not describe 
the confinement-deconfinement transition. 

The Lagrangian density of the chiral quark-meson 
model is given by 

C = \ {d^f + qi$q + gqMq + U(o, tt), (1) 

where the 0(4) representation of the meson fields is 
4> = (c,7?) and the corresponding SU(2)l ® SU(2)r chi- 
ral representation is given by M = o+ir-TVjs. There are 
Nj = 4 mesonic degrees of freedom coupled to N j = 2 
quark flavors q. 

The potential U(o, 7?) is given by 

U(o, tt) = V + ^ - co. (2) 

In vacuum, the 0(4) symmetry of the Lagrangian ([11 is, 
for m 2 < 0, spontaneously broken to 0(3). This leads to 
a nonvanishing scalar condensate [3| (o) = f n . The ex- 
plicit symmetry breaking term co in the potential gives 
the pion a mass m v . At the three level, cq = f-Trrn 2 yields 



the physical pion mass in vacuum. For convenience we in- 
troduce a dimcnsionless parameter h = c/cq, the reduced 
external field, as a measure for the strength of the sym- 
metry breaking term. In a medium, the chiral symme- 
try of the Lagrangian is restored, leading to a vanishing 
scalar condensate at some critical temperature and/or 
density. 

To explore the critical properties near the chiral phase 
transition we use the Functional Rcnormalization Group. 
The FRG is an important tool for addressing nonpertur- 
bative problems within quantum field theory The start- 
ing point is the infrared (IR) regularized effective average 
action Tj. [3, H3, [H, [2{|. This is a generalization of the 
standard effective action, where one-particle modes with 
a momentum less than a scale parameter k are suppressed 
by means of a regulator function. 

The FRG equation for the effective action reads [bHI^ 

d fe r fc [$>] = ^Tv^d k R kB (r { p[^^j]+R kB ^ I 

-Tr (d k R kF (4 2) [$, V] + R kF ) , (3) 

(2) 

where T k denotes the second functional derivative of 
Ta;[ < I ) ,'0] with respect to the field variables and corre- 
sponds to the inverse two-point function at the scale k. 
The trace in Eq. ^ denotes a momentum integration 
and a summation over all internal indices like e.g. flavor, 
color, and/or Dirac indices. The fields $ and ip denote 
bosonic and fcrmionic fields, respectively. 

The effective average action T k governs the dynamics 
of a theory at the momentum scale k and interpolates 
between the bare action T k= \ = S and the full quantum 
effective action r^ = o = T. The regulator function R k , 
which suppresses the small momentum modes, is to some 
extent arbitrary. The derivative of this function dtR k is 
peaked at the scale k and thus implements the Wilsonian 
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idea of successively integrating out shells in momentum 
space. 

The effect of fluctuations is thus gradually included 
in Tk by solving the FRG flow equation. In Euclidean 
space-time 

+ + ■ tt^) ijj + U k (p, cr)} , (4) 



where the field p is given by 



2 V 2 V 



(5) 



and the fermionic field ip carries two flavors correspond- 
ing to the up and down quarks. The matrix 7 is the 
Euclidean analogue of 75. 

In the following we assume that each spectral function 
is dominated by a pole, corresponding to a quasiparti- 
cle. Furthermore, we neglect the wavefunction renor- 
malization for both the bosonic and fermionic fields 
(Z^f. = Z^^ = 1), i-e. the anomalous dimension is 
set to zero. Thus, changes of quasiparticlc properties 
in the medium are accounted for, but the fragmentation 
of singcl-particlc strength is ignored. Finally, we neglect 
the scale dependence of the Yukawa coupling g in Eq. 
(jl]). Consequently, the only scale dependence we are left 
with is that of the potential Uk(p)- 

This approximation to the average effective action 
in ([4]) corresponds to the leading order derivative expan- 
sion, which is obtained as the leading term in a system- 
atic expansion in powers of derivatives of the fields. For 
uniform field configurations, the average effective action 
Tfe evaluated in this approximation is proportional to the 
effective potential U k (p,a) 1 since 



T k = / d 4 xU k (p,a). 



(6) 



The leading order derivative expansion offers a relatively 
transparent framework, which nonetheless yields a very 
useful description of the critical fluctuations near a phase 
transition. 

In addition to defining the form of the effective action, 
one also needs to specify the regulator functions R k - If 
suitably chosen, the momentum integration is both in- 
frared and ultraviolet finite, leading to a numerically sta- 
ble solution of the flow equations. We use the so called 
optimized regulator functions [30T ]. 



for bosons and 

R7?M = i 




- 1 9{k z - q z 



(7) 



(8) 



for fermions. Such regulator functions in general lead to 
a stabilization of the RG flow and an improved conver- 
gence towards the physical theory [3(3 . Due to the theta 



function, the regulator acts as a momentum-dependent 
mass term for all modes with q 2 < k 2 . 

A. Flow at finite temperature and chemical 
potential 

In the previous section we have introduced the basic 
concepts of FRG. In the following we will use this method 
to find the nonperturbative thermodynamic potential at 
finite temperature and density. 

We treat bosons and fermions in thermal equilibrium in 
the standard imaginary time formalism. In thermal equi- 
librium the (fermion) boson fields satisfy (anti-) periodic 
boundary conditions in the Euclidean time direction with 
periodicity 1/T. The momentum integration is replaced 
by a Matsubara sum as follows: 



d d q 
(27r) d 



jd-l 



<1 



{2ir) d - 



(9) 



where 



q (n) = 2mrT, q (n) = (2n + l)nT (10) 



for bosons and fermions, respectively. 

We employ a modified form of the optimized regula- 
tor, where the Euclidean 4-momcntum squared q 2 in Eqs. 

and ([8]) is replaced b y th e 3-momentum squared q 2 . 
Thus, for bosons we use [31| 



<L(q) = (fc 2 -q 2 Mfc 2 -q 2 )- 



nil 



The finite quark chemical potential p is introduced by 
the following substitution of the time derivative 



do -* d + ip 



(12) 



in the fermionic part of the effective action (|4]) . At finite 
p, the fermionic regulator function is then of the form 



r °f1(i) = {i + nn 




(g0 + ^ )2 + fc "-iW-q 2 ). 



(90 + i^) 2 + q 2 



(13) 

Clearly, the modified regulators are invariant only un- 
der spatial rotations of the momentum q in a particular 
frame, but not under Euclidean rotations involving the 
imaginary time direction. However, this is not a cru- 
cial issue in calculations at finite T and p, where the 
heat bath defines a preferred frame. Furthermore, the 
Euclidean invariant form leads to various difficulties, as 
pointed out in [3l| . A great advantage of the regulators 
(|11|) and (|13p is that, in the quasiparticlc approximation, 
the Matsubara sums for the one-loop diagrams are iden- 
tical to those appearing in a free theory. 

Thus, with the regulators (jTTJ) and (fT3|) one obtains the 
following flow equation for the thermodynamic potential 
density (tt = (T/V)T) 
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FIG. 2: Scaling of the order parameter with t and h. The logarithmic fits yield the critical exponents j3 (left panel) and 8 (right 
panel) . 



d k fl k (T, /i) 



12tt 2 



E„ 



l + 2n B (E n ) 



E„ 



where n B {E 7r ^ a ) and n F {E q ), n F {E q ) are the bosonic and 
fermionic distribution functions 



n F (E q ) 
Here 

Ei, 



1 

e (E q -n)/T + 1 - 



n F {E q ) 



1 



e (E q +n)/T + i 



W k , E a 



k 2 + ti> k + 2pQ.'l 



is the pion and sigma energies, where fl' k and 51//' denote 
the first and the second derivatives of the thermodynamic 
potential with the explicit symmetry breaking term re- 
moved, = Qk + pk, with respect to p and 



Eg = V/P + VP 

the quark energy. Finally, the quark degeneracy is given 
by v q = 2N c N f . 

The flow equation (fT4"|) is identical to the one obtained 
within the proper time renormalization group (PTRG) 
scheme [H, EH. In fact, it was shown by Litim [3(| that, 
in vacuum, the PTRG scheme, with a properly chosen 
cutoff function, is equivalent to the optimized FRG flow 
in the leading order derivative expansion. 

The flow equation for Q^iT^p) can be solved numeri- 
cally by either discretizing the potential on a grid 0] or 
expanding it in powers of the fields (tr, p). In this paper, 



l + 2n B (E a ) 



r 



2//, 
~E, 



q -[l-n F {E q )-n F {E q ) 



(14) 



we employ the second method and expand the poten- 
tial flk(T,p) in a Taylor series around the minimum at 
Co = V^PO 



Z — ^ ?TJ 1 



(15) 



The symmetry breaking term is treated explicitly, leading 
to 



fi fc (T) = n fc (T) - ca h . 
The stationarity condition 

an* (T). 



da 



0. 



(16) 



(17) 



determines the position of the physical minimum at the 
scale k. Thus, for c ^ 0, i.e., for an explicitly broken 
symmetry, we find 



c = aicr 



(18) 



which relates the coupling a\ and the expectation value 
of the scalar field a . 

We truncate the expansion in Eq. (|15[) at m = 3. Using 
the identity [H[ 



dn' k _ dn' k dp t dn' k 
~~d~k 



(19) 



dp dk dk 

and the corresponding relations for Cl' k ' and SI'S', we then 
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obtain the flow equations 

d k a = — d k p k + d k Sl k , (20a) 
\J 2pk 

dkPk = —( ho hrz \ d * Q 'k> ( 20b ) 

d k a 2 = a 3 d k p k + d k tt' k \ (20c) 
d k a 3 = d k Q'l', (20d) 

where we have introduced the shorthand notation d k = 
d/dk. In Eq. (fT9)) the terms dO.' k /dp etc. are evaluated 
using the expansion (|15| truncated at m = 3. 

In the following we consider only the case of vanishing 
net baryon density, i.e., p = 0, while the external field c is 
kept finite. Results for the chiral limit are then obtained 
by letting c — ► 0. 

The flow equations are solved numerically starting at 
a cutoff scale A = 1.2 GcV. The initial values for the 
parameters are chosen such that the physical pion mass 
jtJjt = 138 MeV, the pion decay constant f n = 93 MeV 
and the sigma mass m a = 600 MeV are reproduced in 
vacuum, using the regulators (fTTj) and (fl~3|) . for a con- 
stituent quark mass of m q = 300 MeV. The strength of 
the external field Co and the Yukawa coupling (g = 3.2) 
are fixed by the pion mass and the constituent quark 
mass. The pion and the sigma masses, at a given mo- 
mentum scale k, are given by 

m l,k = ~m=> m l,k = -75= + 2p k a 2 , k , (21) 

while the constituent quark mass is obtained from 

m* fc = 2g 2 Pk . (22) 

The numerical solution of the flow equations yields a 
non-perturbative thermodynamic potential at finite tem- 
peratures. This can then be used to explore the critical 
properties of thermodynamical observables in the vicinity 
of the chiral phase transition. 

III. CRITICAL BEHAVIOR AND 0(4) SCALING 

It is well known that the effect of critical fluctuations 
on the properties of a system close to a second-order 
phase transition can be efficiently computed by means of 
the Wilsonian renormalization group [321 ]. A systematic 
account of fluctuations is essential for a realistic descrip- 
tion of the critical properties of such systems. The FRG 
approach offers a powerful tool for computing not only 
universal quantities, like the critical exponents, but also 
the full thermodynamics of a strongly interacting many- 
body systems. Here we focus on the critical exponents of 
the chiral phase transition. 

In the chiral limit, the chiral quark- meson model ex- 
hibits a second order phase transition. The critical be- 
havior of this model (and of QCD) is expected to be gov- 
erned by the fixed point of the 0(4) Heisenberg model 



in three dimensions. Universality class arguments lead 
to predictions for the functional form of various thermo- 
dynamic quantities in the vicinity of the critical temper- 
ature T c , the critical exponents. These all emerge from 
the scaling form of the singular part of the free energy 
density 

F s (t, h) = b- d F s (b yt t, b Vh h). (23) 

The scaling relation (|23]) determines the properties of T s 
under a scale transformation of all lengths by a factor 
b. There are two relevant parameters in F s (t, h) which 
control the critical behavior: the reduced temperature 
t = (T — T c )/T c and the reduced external field h = c/cq. 
The dimcnsionless quantity T s is the free energy density 
in units of . All physical information concerning 

the phase transition can be encoded in the equation of 
state which relates the order parameter, the external field 
and the reduced temperature in the vicinity of the crit- 
ical point. The leading singularities at the critical point 
are then completely controlled by the thermal and the 
magnetic critical exponents, yt and y^. 

In this section we compute critical properties and the 
equation of state of the quark meson model using the 
FRG method and confront our results with those of the 
0(4) universality class. We compare the behavior of var- 
ious physical observables in the critical region with that 
expected from the universal scaling function (|23p and ex- 
tract their critical exponents. 

A. Scaling of the order parameter 

A phase transition involving spontaneous symmetry 
breaking, is signalled by the vanishing of an order param- 
eter, when approached from the broken symmetry phase. 
In the quark meson model, a possible order parameter 
related to chiral symmetry breaking is obtained by dif- 
ferentiating the free energy with respect to the external 
field c. Thus, it follows from Eqs. ([2]) and ([5]) that the or- 
der parameter can be identified with the thermodynamic 
average of the sigma field 

(a) =a = ^2p~ . (24) 

In the FRG approach, the value of the order parameter 
depends on the scale k; the flow of (a) follows from equa- 
tion Eq. (|20b[) . The physical value of the order parameter 
at a given temperature T is then given by y/2 Pk-*Q, where 
Pfc^o is the value obtained at the end of the flow (k — > 0). 

In Fig. [1] we show the resulting dependence of the or- 
der parameter on the reduced temperature t for several 
values of the external field h. In the chiral limit, i.e. for 
h = 0, the order parameter (a) smoothly decreases with 
increasing t and finally vanishes at the critical tempera- 
ture t = 0. Such a behavior is characteristic of a system 
with a second order phase transition. For finite h, the chi- 
ral symmetry of the Lagrangian is explicitly broken and 
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FIG. 3: The ratio BIf a \ as a function of t (left), for h = 6 x 10 13 , and as a function of h for three values of the reduced 
temperature t. 



the second order phase transition turns into a crossover 
transition. Thus, for small but finite h, the order param- 
eter (cr) changes fast in a narrow temperature interval, 
at the so called pseudo-critical temperature. However, in 
this case the order parameter remains non-zero, even at 
very high temperatures, as seen in Fig. [T] and there is no 
true phase transition. 

In the vicinity of the critical point the scaling of the 
order parameter (cr) can be obtained from the universal 
form of the free energy. Indeed, by choosing the scale 
factor b in the Eq. flUj such that b Vh h = 1 or b Vt \t\ = 1 
one finds 



and 



T s = \t\ d "T s {t\t\-\h\t 
respectively. Consequently, using 



(25) 



(26) 



(27) 



where (a) = (cr) / f n is the reduced order parameter, one 
finds the corresponding equation of state for the magne- 
tization 



(cr) = hV S F h (z) 



(28) 



and 



(a) = \t\Pr s (t\t\- L ,h\t\-P°) 



-0S\ 



(29) 

with z = th -1 ' 06 and the universal function Fh(z) is 
obtained from Eq. (|25|) by differentiation with respect to 
h. The prime on T a in (|29p denotes the derivative with 
respect to the second argument. 

The well known scaling behavior of the magnetization 



B{-t) ( \ h=0, t<0 
B c h x / S , t = 0, h>0 ' 



(30) 



follows from Eqs. ([281) and p9|l- The constants B c = 
Fh(0) and B = ^(—1,0) characterize the magnetization 
on the so called coexistence line [33| (t < 0, h = 0) 
and at the pseudo-critical point (t = 0, h > 0). In the 
derivation of Eqs. ([25 ]) -([30 l) we have used the definitions 
of the thermal and magnetic critical exponents 



1 1 05 

Vt = -i Vh = — = — , 

V V r V 



as well as the scaling relations 

7 = 1/(2-77), du = (3(1 + 6) 



2 

0' 



and 



a = 2 — dv. 



(31) 



(32) 



(33) 



In the chiral limit, the scaling behavior of the order pa- 
rameter is controlled by the critical exponent f3, when 
the critical temperature t = is approached from below. 
Thus, (3 is given by the slope in a double-logarithmic plot 
of (cr) versus t. The critical exponent 6, on the other 
hand, characterizes the scaling of (cr) as h — > for t = 0. 

In Fig. [5] we show the scaling properties of (cr) close to 
the critical point (t — > 0) for vanishing external field as 
well as the scaling with h at the (pseudo)critical point 
(t = 0). In both cases the order parameter shows the 
expected scaling behavior. This demonstrates that, in 
the parameter range considered, the free energy of the 
quark meson model is dominated by the singular part 
Eq. (f2"3"]) . The critical exponents extracted from a fit to 
the results of Fig. [2] arc 



/3 ~ 0.402 and 



4.818. 



(34) 



Both values are consistent with 0(4) universality. 
Indeed, in a Monte Carlo simulation of the three- 
dimensional 0(4) spin model, Kanaya and Kaya [26| 
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find (3 = 0.3836(46) and 5 = 4.851(22). The agree- 
ment with this presicion calculation is, given our simple 
Ansatz, quite satisfactory. Consequently, we conclude 
that the critical properties of the chiral order parame- 
ter in the quark-meson model are indeed governed by 
the 0(4) universality class. Our results, obtained in the 
relatively simple leading order derivative expansion, illus- 
trate the efficiency of the FRG approach to correctly ac- 
count for a non-perturbativc, long distance physics near 
a phase transition. 

In the mean field approximation one also finds scaling 
of the order parameter (a) close to the critical point, 
but with different critical exponents, /3 = 0.5 and S = 
3. Thus, the singularity at the chiral phase transition is 
substantially modified by fluctuations. 

In the formulation of the flow equation for the quark- 
meson model we have neglected the wavefunction renor- 
malization for both bosonic and fermionic fields. The 
good agreement with the Q(4) lattice results as well as 
with previous FRG studies [19( , where the anomalous di- 
mension was taken into account, indicates that the influ- 
ence of the wavefunction renormalization on the critical 
properties at the chiral phase transition is quite small. 

So far we have discussed the scaling of the order pa- 
rameter (a) with t for h = in the chirally broken phase 
and with h at t = 0. In the chiral limit, the order pa- 
rameter vanishes identically in the symmetric phase, i.e. 
for t > 0. Thus, the order parameter does not exhibit 
a singularity when the critical point is approached from 
above. However, since for h ^ the order parameter is 
always non-zero, it is interesting to explore its behavior 
for small but finite values of h, when the critical point 
is approached from the broken (t < 0) as well as from 
the symmetric (t > 0) phase. The scaling behavior of 
the order parameter is in this case characterized by a so- 
called effective critical exponent (24[. For a given h, this 
exponent is defined by 



d\og (g) 
dlogt 



t d(a) 



(35) 



tent with the value of the critical exponent /3, extracted 
from the scaling properties (|30|) of the order parameter. 

In the right panel of Fig. [3] we show the dependence 
of the effective critical exponent upon h in the broken 
phase for three different values of t. For t < 0, the limit 
h — ► corresponds to z — > — oo and consequently R^ 
again converges to the critical exponent /3. On the other 
hand, the limit h — > oo corresponds to z — > 0. In this 
case, the order parameter (a) ~ h 1 / 8 is, to leading order, 
independent of temperature and the effective critical ex- 
ponent approaches zero. The transition between the two 
regimes takes place at \z\ ~ 1. 

Also in the symmetric phase (t > 0), the exponent R( a \ 
shows a power law behavior, however with a different 
value of the exponent. In the transition region, R^ a \ 
behaves like a smoothened step function, interpolating 
between the asymptotic regions. The transition region 
extends from z ~ -1 to z ~ 1. Hence the width of the 
step in t is expected to scale with ft, 1 /' 35 . The effective 
critical exponent changes from R^ ~ 0.4 in the broken 

phase to i?^ ~ —1.5 in the symmetric phase. This 
change in the critical behavior of the order parameter 
can be understood using the universal scaling behavior 
of the singular part of the free energy. To sec this, we first 
write the equation of state (|28|) and |29|) in the Widom- 
Griffiths form 



y = f(x), 



where 



V 



(36) 



(37) 



The scaling relations ([3D)) are, modulo constant fac- 
tors, obtained if /(0) = 1 and /(— 1) = 0. The constant 
factors are recovered by a simple rescaling of t and h. A 
comparison of the Widom-Griffiths form of the equation 
of state (|3"o) with that given in Eq. (f2"5)) . shows that the 
following relations hold: 



F h (z) = f(x)-^ s 



From Eqs. ([30]) and J35]) it is clear that for h = and Consequently, Eq. (j2"5|) can be rewritten as 



t < 0, the ratio R(~\° coincides with the critical exponent 
(3. For h 0, Eq. (|55|) can be used to extract the lead- 
ing singularity of the order parameter in the symmetric 
region close to the critical point. 

In the left panel of Fig. [3] we show the resulting de- 
pendence of the effective critical exponent (|35|) on the 
reduced temperature for a fixed, very small, value of the 
external field h. The exponent R^ exhibits the expected 
scaling behavior. In the broken (i < 0) as well as in the 
symmetric (t > 0) phase, it is essentially independent of 
t. The effective critical exponents are different in the bro- 
ken and symmetric phases. A weak external field h ~ 
and t < corresponds toz-> — oo in Eq. ([28]). It follows 
that the asymptotic form of the scaling function is given 
by [H| Fh(z) ~ (— z) 13 for z — > — oo. Thus, in the broken 
phase (t < 0) the exponent R 1 ? is, as expected, consis- 



(a) = BchV'fixY 



1/5 



(38) 



(39) 



where we have reinstated the constant factor B c . This 
equation of state can be used to explore the scaling prop- 
erties of (a) when approaching the critical point from the 
symmetric phase. For t > and for h — > 0, the chiral 
order parameter (normalized to its value in vacuum) is 
very small. Consequently, the argument x of the func- 
tion f{x) in equation (pj9"]) is very large and the scaling 
behavior of (a) for t > and for small but finite h is de- 
termined by the asymptotic form of the scaling function 
f(x) for x — » oo. In this limit, f(x) is given by Griffiths' 
analyticity condition [35[, 



oo 

E 



7-2(n-l)/9 



(40) 
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FIG. 4: The susceptibility % a as a function of the reduced 
temperature t for several values of the field h. 



of a susceptibility, e.g. at the critical end point, signals 
a zero-mass mode of the corresponding effective field. 

In the chiral quark-meson model, we are not only deal- 
ing with fluctuations of the order parameter, i.e. of the 
sigma field, but also with a pseudo-Goldstone mode, cor- 
responding to fluctuations of the pion field. The sigma 
and the pion represent the longitudinal and transverse 
modes of the 0(4) field <fi, respectively. One therefore 
distinguishes between the longitudinal (er) and transverse 
(it) order parameter susceptibilities. 

The longitudinal susceptibility is given by the deriva- 
tive of the order parameter with respect to the external 
field 



Xl = Xa 



dh 



(42) 



whereas the transverse susceptibility, which is directly 
related to the fluctuation of the Goldstone modes, is ob- 
tained from 



Retaining only the leading term, f{x) ~ x 7 , one finds 
that for t > and h — ■» the equation of state scales as 



Xt = Xtt 



h 



(43) 



(a) - t">h. 



(41) 



Thus, for t > and h ^ 0, we expect R 1 ?, = —7, in good 
agreement with the results shown in Fig. [3] This result 
implies [34| that for z — > 00, Fi l (z) ~ z -7 . We note that 
in the chiral limit (h = 0), the order parameter vanishes 
identically for t > 0, as expected. 

From the scaling behavior of (a) shown in the left panel 
of Fig. [3j we find 7 ~ 1.6. This is consistent with the lat- 
tice result for the 0(4) spin system obtained in ref. [26T |. 
7 ~ 1.5. Thus, the FRG method applied to the quark 
meson model yields a scaling behavior of the order pa- 
rameter, which is consistent with that found for the 0(4) 
universality class. This is the case at the critical point, 
on the coexistence line and when approaching the critical 
point from the high temperature phase for a finite value 
of the external magnetic field. 



B. Scaling of the chiral susceptibilities 

The critical properties of a thermodynamic system can 
be explored by studying the fluctuations of various ob- 
servables. In particular, the fluctuations of the order 
parameter probe the order of the phase transition and 
the position of a possible critical end point. In statistical 
physics, fluctuations are reflected in the corresponding 
susceptibilities x- I n the case of the chiral order param- 
eter, the corresponding susceptibility is defined as the 
response to a change of the external field h; the suscepti- 
bility x is obtained by taking the second derivative of the 
effective potential with respect to h. The susceptibility in 
other channels is computed analogously. It follows that 
the susceptibilities are inversely proportional to the cor- 
responding mass squared. Consequently, the divergence 



The latter expression is a direct consequence of the 0(4) 
invariance of the free energy for h = and follows from 
the corresponding Ward identity [36| . 

The critical behavior of the longitudinal and trans- 
verse chiral susceptibility can be obtained from the scal- 
ing form of the singular part of the free energy or directly 
from the magnetic equation of state (jST))) . 

In Fig. [4] we show the longitudinal susceptibility Xa as 
a function of the reduced temperature for several values 
of the symmetry breaking field h. At finite h, this suscep- 
tibility shows a peak structure with a maximum at the 
pseudo-critical temperature t maJC (h). With decreasing h, 
there is a shift of t max and an increase of the fluctua- 
tions at the pseudo-critical point; in the chiral limit x<j 
diverges at the transition temperature, because the mass 
of the sigma field vanishes. 

Using the results shown in the Fig. one can define a 
pseudo-critical line in the (t, /i)-plane that characterizes 
the dependence of the transition temperature t max on the 
symmetry breaking field h. Renormalization group ar- 
guments imply that the pseudo-critical temperature and 
the maximum of the longitudinal susceptibility should 
exhibit universal scaling according to [371 ] : 



tu 



(h) ~ T h 1/{l+l3) 



c P tZL, ( 44 ) 



with 0(4) critical exponents. 

In Fig. [5] we confront the scaling behavior of the sus- 
ceptibility with Eq. (|44p . The critical exponents are 
extracted by performing a linear logarithmic fit to the 
pseudo-critical temperature and the maximum value of 
the susceptibility. We find 1/(7 + (3) = 0.494 and 
7 = 1.618, which implies (3 = 0.406, in good agreement 
with the results obtained above using Eqs. (J30|) and (|4Tj) . 
Thus, the scaling of of Xo with the symmetry breaking 
field h is also consistent with 0(4) universality. 
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FIG. 5: Scaling of the pseudo-critical temperature defined by the maximum of the susceptibility \ a (left) and of the maximum 
value of x& (right). 



In the symmetric phase, the longitudinal and the 
transversal susceptibilities coincide (x<r = X-n-)- Conse- 
quently, the strength of the singularities at the critical 
point t = are controlled by the same critical exponents. 

To find the scaling behavior of \ when approaching 
the critical point from the symmetric phase we use the 
Widom-Griffiths form of the magnetic equation of state 
([55)) together with the large-x expansion of the scaling 
function f{x). 

From Eqs. (gTJ), (gj) and (O it is clear that for h -> 
in the symmetric phase (t > 0) 

X a = X, ^ t~\ (45) 

Thus, when approaching the critical point from the sym- 
metric phase, the scaling behavior of the longitudinal and 




-12 -11 -10 -9 

log(t) 

FIG. 6: Double logarithmic plot of the longitudinal suscepti- 
bility in the symmetric phase for a small h as a function the 
reduced temperature. 



transverse susceptibility is controlled by the critical ex- 
ponent 7. 

In Fig. [6] we show the scaling behavior of the suscep- 
tibility for t > and in the limit of h — > 0. From the 
slope of the line we obtain 7 = 1.575, in good agree- 
ment the value obtained from the analysis of the scaling 
of the chiral susceptibility Xa following Eq. ([44]) . The 
small difference is may be due to the uncertainty in the 
determination of the critical temperature. 

In the broken phase, on the coexistence line, the crit- 
ical behavior of the longitudinal and transverse suscep- 
tibilities differ. In the chiral limit, the order parameter 
is finite in the broken phase as long as t < 0. Thus, 
for T < T c the transverse susceptibility (|4"5|) diverges in 
the chiral limit, x-k ~ due to the appearance of 

Goldstone bosons. A less obvious result is the divergence 
of the longitudinal susceptibility on the coexistence line, 
where to leading order Xa ~ h~ 1 / 2 0, [HI]. 

The scaling properties of the susceptibilities can be 
obtained directly from the Widom-Griffiths form of the 
magnetic equation of state. Indeed, using Eqs. (|59"]) . (|4"1?|) 
and (|43j) one finds 

y -Bk l/S - 1 IMiEt^H (46) 

for the longitudinal and 

X, = B^/^fix)^ 5 , (47) 

for the transverse susceptibility. 

In the vicinity of the coexistence line, where x — > — 1, 
the scaling function is, in a three-dimensional system, 
found to be of the form [H, [38| 

f(x)^c f (l + x) 2 . (48) 

Using (gSJ) in (gBJ) and ([S7)l. one then finds the critical 
behavior of the susceptibilities for t < and h — > 0: 

Xa - ft-V2(_f)/3-W/2 (49) 
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FIG. 7: Scaling of Xo- and A> in the broken phase (t < 0) as functions of h (left) and t (right) 



and 



(50) 



The divergence of the susceptibilities on the coexistence 
line is due to the transverse modes (33j , whose mass term 
is proportional to the symmetry-breaking term, m^. ~ 
h. The different exponents for x<r an d X-n reflect the 
direct and indirect dependence of the susceptibilities on 
the soft transverse modes. The critical behavior of the 
longitudinal susceptibility close to the coexistence line, 
fixes the exponent of the sublcading term of the scaling 
function defined in (l28l) for z < 



F h (z)~(-zf(l + C(-z)-W) 



(3S/2-, 



(51) 



where C is a constant. 

In Fig.[7|we show log-log plots of the longitudinal and 
transverse susceptibilities, obtained by solving the flow 
equations for the quark meson model (|20j) . as functions 
of h and t near the coexistence line. A linear fit to the 
results shown on the left yields 0.499 and 1.000 for the 
critical exponents of x<y an d Xv> m excellent agreement 
with the expected scaling behavior with h implied by Eqs. 
(|4"9")l and ([50")) . The corresponding fit to the i-dependence, 
shown in the right plot, gives (f3 — /3S/2) = —0.617 and 
/3 = 0.395, which implies 8 = 5.124, in fair agreement 
with our previous results for the critical exponents (3 and 
8 obtained from the scaling of the order parameter. 

In order to exhibit the scaling behavior of the suscepti- 
bilities in a broad parameter range, which covers both the 
broken and symmetric phases, we introduce the ratios 



Xt, 



dt 



(52) 



for the transverse x-k and the longitudinal Xa susceptibil- 
ities. The ratios i2£ define effective critical exponents, 
analogous to Eq. (f35| . 



In Fig. [8] we summarize the t and h dependence of the 
and in the vicinity to the critical point. On the 
left we show how the critical exponents of the suscep- 
tibilities interpolate between the broken and symmetric 
phases, for fixed h, while on the right the ^.-dependence 
of i?£ a in the symmetric phase is illustrated for a few 
values of t. 

In the symmetric phase, x-n arid x<j coincide and their 
critical behavior is determined by the same exponent 
(7 ~ 1.58), as seen in Fig. [5]Tcft. This corresponds to 
the limit z — > 00, and is consistent with the effective 
critical exponent obtained for the order parameter for 
t > 0. The dependence upon h in the symmetric phase, 
shown on the right, is also identical for the longitudinal 
and transverse susceptibilities. Here the z — ► 00 limit 
corresponds to h — > for fixed t > 0. Thus, also in 
this limit, R^ ^ — > —7, as shown in Fig. [5]-right. Again, 
the h — > 00 limit corresponds to z — > 0, where the or- 
der parameter and consequently also the susceptibility 
approaches a temperature independent function of h and 
the corresponding effective critical exponent converges to 
zero with increasing h. The transition between the two 
regimes again takes place at z ~ 1. 

In the broken phase (see Fig.[8]-left), the critical expo- 
nents differ, R^ ~ 0(1 — 8/2) and R^ ~ /3, in complete 
agreement with Eqs. (|4"9")l and (|50"|) . 

In section fill Al we determined the critical exponent 8 
by studying the scaling of the order parameter (a) with 
h for t = 0. We note that this critical exponent may 
be determined in an alternative way, by exploiting the 
scaling properties of the susceptibilities [25|, [39|, Ho| X™ 
and Xa- Consider the ratio 



A(t,h) 



X.7T 



h d(a) 
W)~dh 



(53) 



which has the useful property that its value at t = is 
independent of h within the scaling region. Indeed, using 
the parametrization of the susceptibilities (J3U) and ([47)1 



11 




as well as the properties of the Widom-Griffiths scaling 
function f(x) discussed above, one finds 

f 1, t>0, 
lim A(t,h) = { 1/6, t = 0, (54) 
{ 0, t<0. 

Consequently, the ratio A(t, h) can be used to directly 
extract the value of the critical exponent 6, without re- 
sorting to logarithmic fits. Moreover, this ratio simulta- 
neously reveals the position of the critical point, since in 
the scaling regime, all A(t,h), plotted as functions of t 
for different h, cross at t = 0. 

In Fig.[S]-left we show A(i, h) as a function oft for four 
values of h. As expected all lines cross at t = 0, yielding 
the value of the critical exponent 6 « 4.91. This value is 
in good agreement with the one obtained above, using a 
logarithmic fit to the dependence of the order parameter 
on h. 

In the right panel of Fig. [5] we show the scaling of A as 
a function of h, for fixed t > 0. In the chiral limit, that is 
for h — 0, the longitudinal and transverse susceptibilities 
coincide above the critical temperature. Consequently 
A = 1 in the limit h — > 0, as shown in the plot. In the 
opposite limit, i.e. for h — ► oo, we find A = 1/8. The 
coincidence of the two limits h — > oo and t — ► follows 
from the explicit form for this ratio in terms of the scaling 
equation of state ((2SJ), A = 1/6- {z / j35){F' h (z) / F h {z)). 
Since the limits t — > for finite h and h — > oo for finite 
t both imply z = th.- 1 /' 35 ->■ 0, the result A = is 
obtained in both cases. 



C. Correlation lengths and their critical behavior 

A correlation length £ governs the exponential decay 
of the corresponding correlation function with increasing 
separation. The critical behavior of correlation lengths 



is an important characteristic of phase transitions. At 
the critical point of a second-order phase transition, the 
order-parameter correlation length diverges, heralding 
the onset of a long-range order. 

The decay of a correlation function is controlled by the 
relevant mass scale in the system. In a chiral model, such 
as the one considered here, there are two such scales, the 
transverse mass connected with the pion and the longitu- 
dinal mass related to the sigma field. Consequently, like 
in the case of magnetization, one introduces a transverse 
and a longitudinal correlation length £„■ and £ a . In the 
critical region for t > 0, the correlation lengths are di- 
rectly related to the corresponding susceptibilities @, El| 

X*~£ ; 2 , Xt~e t - (55) 

Consequently, in the chiral quark-meson model, the cor- 
relation lengths are determined by the pion and sigma 
masses 

, £r~— ■ (56) 

In the broken symmetry phase, the soft transverse 
modes, which renormalize the longitudinal susceptibility 
in a non-trivial way, as shown in Eq. 1491 also profoundly 
influence the longitudinal correlation length 42]. This 
modifies the relation between the correlation length and 
the susceptibility in the longitudinal channel for t < 0. 
We do not address this problem here. 

In Fig. [TO] we show the correlation lengths, £ CT and 
at the critical point t = as functions of the external 
field h. As expected both the longitudinal and transverse 
correlation lengths diverge as h — > 0. The divergence of 
£ CT reflects the softening of order parameter fluctuations 
as the critical point is approached, while that of ^ signals 
the appearance of Goldstone bosons at the second order 
phase transition. 

The scaling of the correlation length in the critical re- 
gion is obtained from the scaling function (|28|) and the 
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FIG. 10: The correlation lengths £ CT and f w at the critical 
point t — as functions of the external field h. 



relation ((55)) . On the critical line (t = 0, h — ► 0), the 
transverse and longitudinal correlation length exhibit the 
same singularity, controlled by the exponent v c Q 



(57) 



Their ratio, ^n/^a, is at i = approximately equal to 2 
(see Fig. HU]) . Thus, the value expected for this ratio in 
the broken symmetry phase [421 is obtained also on the 
critical line, in agreement with 7j. 

In the symmetric phase (t > 0), the correlation lengths 
£cr and £ w coincide and the scaling with t is governed by 
the critical exponent v 



£cr,7T ^ t 



(58) 



In Fig.[TT]we show the scaling behavior of the correlation 
lengths on the critical line and in the symmetric phase. 
A fit of the critical exponents yields v c = 0.396 and v — 
0.787, respectively. 



Our results for the various critical exponents are sum- 
marized in Table |U We give the values obtained in the 
logarithmic fits shown in Figs. l2l6 landfTTl which we deem 
to be more accurate than the other determinations. The 
value for the coefficient a is obtained by using the hy- 
perscaling relation (|3"3")) . Our results are compared with 
recent lattice results and mean-field theory exponents for 
the 0(4) model. The mean-field values clearly deviate 
from the FRG and lattice results, showing the impor- 
tance of fluctuations near the critical point. Our ex- 
ponents, obtained with the FRG approach, are in good 
agreement with the Monte Carlo lattice results [26| in 
spite of the fact, that we have neglected the anomalous 
dimension in the flow equations. 

The finding that the anomalous dimension is fairly 
unimportant for the critical behavior, is consistent with 
the relatively small value of the corresponding critical ex- 
ponent r\ in the 0(4) universality class. In spite of the 
fact that we neglect momentum dependent couplings in 
the RG flow, a fairly reasonable value for the critical ex- 
ponent ry, can be extracted by means of the hyperscaling 
relation 



d+2 



V 



d-2 + -q 



(59) 



Using our value for S, we find rj w 0.031, which should 
be compared with the lattice result, rj = 0.0254 [2(|. We 
note, however, that using the middle scaling relation in 
(|3"2")l to determine 77 fails; the resulting value is small and 
negative. This illustrates the uncertainty introduced by 
the leading order derivative expansion. 



IV. SUMMARY AND CONCLUSIONS 

We have discussed the critical properties of the chi- 
ral quark-meson model at the finite-temperature phase 
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FIG. 11: The correlation lengths £ CT and ^ as functions of h at t — (left) and as functions of t at h = (right). As discussed 
in the text, the scaling of the correlation lengths is governed by the critical exponents v c and v. Note that the normalization 
of the correlation lengths differ from that used in Fig. 1101 



TABLE I: The critical exponents of the 0(4) model obtained in lattice calculations, in the mean-field approximation and in 
the present FRG approach. 





P 


S 


7 


V 




a 


0(4) lattice [26] 


0.3836 


4.851 


1.477 


0.7479 


0.4019" 


-0.244 


Mean Field 


0.5 


3 


1 


0.5 


1/3 





FRG (our work) 


0.402 


4.818 


1.575 


0.787 


0.396 





"Obtained using the definition of v c (J3TJ . 
''Obtained using the scaling relation (133 P 



transition. The focus was on the universal scaling behav- 
ior of various physical quantities. Fluctuations and non- 
perturbative effects were accounted for by employing the 
Functional Renormalization Group (FRG) method. We 
derived the flow equation for the scale dependent ther- 
modynamic potential at finite temperature and chemical 
potential in the presence of an external field. The flow 
equation was then solved numerically, by employing a 
truncated polynomial expansion around the minimum of 
the effective action. We have shown that the choice of 
regulator functions (fTTj) and (fl"3"]) leads to transparent ex- 
pressions and reliable numerical results. In this work, we 
have neglected the bosonic and fermionic wave function 
renormalization and the running of the Yukawa coupling. 
We argue that, in the critical region, these effects lead to 
sub-leading corrections to the scaling behavior. 

We have computed various thermodynamical quanti- 
ties that exhibit critical behavior near the chiral phase 
transition. In particular, we have analyzed the scaling 
behavior of the chiral order parameter, its longitudinal 
and transverse susceptibilities as well as the correlation 
lengths in the chiral limit and in the presence of an ex- 
ternal symmetry breaking field. We found that the chi- 
ral order parameter and its susceptibilities scale near the 
chiral phase transition following the universal scaling be- 



havior of the singular part of the free energy. Here we 
have extracted the corresponding critical exponents at 
the critical point and at the coexistence line and have 
also analyzed the scaling of the pseudocritical tempera- 
ture and the maxima of the longitudinal susceptibility in 
an external field. 

We have also analyzed the effective critical exponents, 
which govern the leading scaling behavior above and be- 
low the critical temperature for a system in an external 
symmetry breaking field. In this context, the longitu- 
dinal and transverse susceptibilities are of particular in- 
terest, since they diverge at the coexistence line. The 
effective critical exponents obtained numerically, are in 
full agreement with the Widom-Griffiths form of the mag- 
netic equation of state. 

The critical exponents obtained in this work within the 
quark meson model by employing the FRG method, are 
in a very good agreement with the recent lattice results 
for the 0(4) spin system. This confirms the expecta- 
tion that the chiral phase transition of the 2-flavor effec- 
tive chiral Lagrangian belongs to the 0(4) universality 
class. Furthermore, our results demonstrate the power 
of the FRG approach and the validity of the truncation 
approximation to the flow equation. This approach of- 
fers an efficient framework for describing fluctuations and 
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non-pcrturbativc long-distance dynamics near the chiral 
phase transition. 

We note in closing, that the critical region (the region 
of true non-trivial critical behavior) of the chiral quark- 
meson model turns out to be very small in our calcula- 
tion. This is reflected in the very small values for t and 
h needed in order to obtain stable results in the logarith- 
mic fits of the critical exponents. Our finding is consis- 
tent with the suppression of the Ginsburg region [43| in 
O(N) models for large N. However, because the size of 
the critical region is not universal, the quantitative im- 
plications for the chiral transition in QCD are uncertain. 
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